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ABSTRACT 

The fundamental plane (FP) is a widely used tool to investigate the properties of 
early-type galaxies, and the tight relation between its parameters has spawned several 
cosmological applications, including its use as a distance indicator for peculiar veloc¬ 
ity surveys and as a means to suppress intrinsic noise in cosmic size magnification 
measurements. Systematic trends with the large-scale structure across the FP could 
cause serious biases for these cosmological probes, but may also yield new insights 
into the early-type population. Here we report the first detection of spatial correla¬ 
tions among offsets in galaxy size from an FP that explicitly accounts for redshift 
trends, using a sample of about 95,000 elliptical galaxies from the Sloan Digital Sky 
Survey. We show that these offsets correlate with the density field out to at least 
10h -1 Mpc at 4 a significance in a way that cannot be explained by systematic errors 
in galaxy size estimates. We propose a physical explanation for the correlations by 
dividing the sample into central, satellite, and field galaxies, identifying trends for 
each galaxy type separately. Central (satellite) galaxies lie on average above (below) 
the FP, which we argue could be due to a higher (lower) than average mass-to-light 
ratio. We fit a simple model to the correlations of FP residuals and use it to predict 
the impact on peculiar velocity power spectra, finding a contamination larger than 
10% for k > 0.04/i/Mpc. Moreover, cosmic magnification measurements based on an 
FP could be severely contaminated over a wide range of scales by the intrinsic FP 
correlations. 

Key words: methods: data analysis - cosmology: observations - galaxies: distances 
and redshifts - galaxies: fundamental parameters - gravitational lensing: weak - large- 
scale structure of Universe 


1 INTRODUCTION 


Large-scale structure cosmology has moved beyond the ba¬ 
sic measurement of galaxy positions used for example to 
infer galaxy clustering and baryon acoustic oscillations. 
There is much interest in applications that rely on measure¬ 
ments of morphological properties of galaxies. An example 
that has reached a relatively ma ture state is weak gravi¬ 
tatio n al lensing (for reviews, se e iBartelmann fe Schnei der! 

“3T 


l200ll ; iMassev, Kitching fc Richardl l20ld : Kilbingerl 2014lf . 
which typically uses galaxy ellipticities to infer tiny but co¬ 
herent shape distortions, or shear. Analogously, the sizes 
of galaxies can be used to measure weak lensing magni- 
ficati on, which is complementary to shear measurement s 
(e.g., ISchmidt et al.1120121 : iHeavens, Alsing fe Jaffel l2013ll . 
but suffers from a larger shot noise contribution due to 
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the broad intrinsic size distribution (lAlsing et al.l l2014l l. 
This noise term can be effectively suppressed by taking ad¬ 
vantage of the existence of a well-defined relationship be¬ 
tween galaxy siz es and other observables such as the fun¬ 
damental plane dBertin fe Lombardill200d ). The fundamen¬ 
tal plane (FP hereafter) corresponds to a tight relation be¬ 
tween the size, velocit y dispersion, and surface brightness 
of early-type galaxies dDiorgovski fe Pavia 1 19871 1 and has 
been widely used to study the properties of ellipticals. In 
the context of weak lensing magnification the usefulness 
of the FP wa s recen tly demonstrated observationally by 
iHuff fe Graved (l2014li . using a photometric analogue. 

Mapping out the peculiar velocity field (e.g., 
ISpringob et al.l i2014l l is another example of the ex¬ 
ploitation of the FP for the pur poses of a cosmological 
probe (see IStrauss fe Willickl Il995l for a review). Both ap¬ 
plications involve measuring the observed galaxy size, and 
predicting a galaxy size based on the FP. The comparison 
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between the prediction and observation is used to measure 
the quantity of interest, either the size change due to 
lensing magnification, or the line-of-sight peculiar velocity, 
which modifies the redshift and thus the angular diameter 
distance used to translate an angular size to a physical size. 

Both of these cosmological probes make the important 
fundamental assumption that deviations in galaxy sizes from 
the FP relation do not have some underlying correlation 
with the large-scale distribution of matter, other than the 
one induced by the cosmological signal. Such intrinsic size 
correlations are not completely unreasonable though, consid¬ 
ering that a de nsity dependence of FP parameters has been 
observed (e.g. lLa Barbera et ahl [20 101. and that intrinsic 
galaxy cllipticity correlations with the cosmic density fi eld 
are a well-established physical phenomenon (iTroxel fc Ishakl 
l2014ll . Correlations of FP residuals with the matter density 
mimic the correlations between peculiar velocities or mag¬ 
nified galaxy sizes derived from the FP, so that the intrinsic 
correlations in the FP can bias the cosmolo gical inference 
if not accounted for. lStrauss fc Willickl (1 19951) discussed the 
assumption of the universality of distance indicator relations 
used for peculiar velocity studies, including the FP, but con¬ 
cluded that, at the time, no compelling evidence for spatial 
or environmenta l variations existe d. 

Recently. ISaulder et al.l (I2013I . S13 hereafter) used data 
from the Sloan Digital Sky Survey (SDSS) to determine an 
FP for nearly 10 s early-type galaxies, with a plan of en¬ 
abling measurements of the peculiar velocity field. We use 
this sample of unprecedented size to investigate the existence 
of spatial correlations of the FP residuals in the galaxy size 
parameter, and quantify the extent to which such correla¬ 
tions contaminate peculiar velocity and cosmic magnifica¬ 
tion analyses. 

This paper is structured as follows: after describing in 
Section[2]the dataset that we are going to analyse, we present 
FP fits in Section[3] with particular focus on redshift trends. 
In Section^ we summarise the estimators and error mea¬ 
surements for the two-point statistics we employ. Section[5] 
contains our correlation measurements, their modelling and 
interpretation, as well as a discussion of systematics tests. 
In Section[6]we predict the impact of the FP residual corre¬ 
lations on peculiar velocity measurements and cosmic mag¬ 
nification, before concluding in Section[3 

Throughout we assume a spatially flat ACDM cosmol¬ 
ogy with = 0.27, fib = 0.05, crs = 0.8, n s = 0.96, and 
Ho = 70km/s/Mpc unless quantities are specified in terms 
of h with Ho = h lOOkm/s/Mpc. 


2 DATA 


All data used in this paper originates from the Sloan 
Digital Sky Survey (SDSS). The SDSS liYork et al.l [ioool'i 
imaged roughly n steradians of the sky, and followed 
up approximately one million of the detected objects 
spectroscopically, with the relevant spectroscopic sam¬ 
ple for this paper be ing t he flux-limited Main galaxy 
sample (IStrauss et al.l [2002). The imaging was carried 
out by drift-scanning the sky in p hotometric conditions 
(Hogg ct al. 2001; Ivczic ct al. 2004), n five bands ( ugriz ) 
( Fukugita et al.l 19961 : Smith et al .1120021) using a specially- 
designed wide-field camera ( Gunn et al.l I1998h . All of the 
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Figure 1. Distributions of galaxy samples in the rest-frame M r 
magnitude - redshift plane. The heat map (linear scale) corre¬ 
sponds to the full S13 sample, with the black solid contour indi¬ 
cating the densest region containing 68 % of the galaxies. Corre¬ 
sponding contours are shown for the VAGC sample (green), the 
S13 brightest group galaxy (BGG) sample (red), the S13 satel¬ 
lite sample (yellow), and the S13 sample of field galaxies (blue). 
The vertical black dashed line indicates the minimum redshift cut 
imposed for the analysis of two-point statistics. 


data were processed by completely automated pipelines 
that detect and measure photometric p roperties of objects , 
and astrometrically calibrate the da ta ([Luoton et al.ll200ll : 
IPier et aLll2003l : iTucker et alj|2006l ). The SDSS-I/II imag- 
ing surveys were com pleted with a seventh data release 
llAbazaiian et all 120091 ). though this work will rely as well 
on an improved data reduction pipeline that was part of 
the_eighth data release (DR8), from SDSS-III (jAihara et alj 
im3); and an im proved p hoto metr ic calibration (‘ubercali- 
bration’, IPadmanabhan et al.ll2008i ). 

We reproduce the galaxy sample selection of S13 from 
SDSS DR8 as well as their derived galaxy properties. 
In short, we include galaxies with spectroscopic redshifts 
which h ave b een classified as elliptical by Galaxy Zoo 
(jLintott et al.|[201ll'l with high confidence, are better fitted 
by a de Vaucouleurs profile than an exponential profile in 
all five passbands, and have an axis ratio > 0.3. We also 
follow S13 in determining rest-frame absolute magnitudes 
from the de Vaucouleurs model fits, a circularised physical 
galaxy radius, Ro, based on the de Vaucouleurs half-light ra¬ 
dius, an estimate of velocity dispersion, <to, corrected for the 
effects of fixed fibre diameter, and a dimming-corrected sur¬ 
face brightness, /ro, expressed in terms ofi log Jo = —0.4/ro- 
After removing a small number of duplicates and obvious 
outliers in the photometry, and following the trimming of 
the red sequence proposed by S13, we are left with about 


1 We use log to denote the logarithm to base 10, and In for the 
natural logarithm. 
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95, 000 galaxies in the range 0.01 < z < 0.2 usable for FP 
measurements (referred to as the S13 sample henceforth). 
We focus in our analysis on the r and i bands as these have 
the best photometry as well as the tightest FP relations. 

M oreover, we also use t he NYU value-added galaxy cat¬ 
alog jBlanton et al.l 1200-51 ) for SDSS DR7 to construct a 
sample that will serve as a tracer of the matter distribu¬ 
tion (559,028 galaxies, referred to as the VAGC sample). We 
also match the galaxies in t he S13 sample with the updated 
SDSS DR7 group catalog llYang et al.ll2007l . 120121 1 to iden¬ 
tify the satellites, BGGs (brightest group galaxy, treated as 
central galaxy) and field galaxies (galaxies in groups of mul¬ 
tiplicity one). We obtain 18944 (19.95%) satellites, 18416 
(19.4%) BGGs, and 54762 (57.7%) field galaxies in the S13 
sample. For 2812 (2.95%) galaxies we did not find a unique 
match in the group catalog, so those galaxies are not rep¬ 
resented in figures that divide the sample based on envi¬ 
ronment. Figure[l] shows the joint redshift and rest-frame 
magnitude distributions of the samples used in this work, 
including the environment subsamples defined above. 


3 FUNDAMENTAL PLANE 

We employ a weighted linear least squares algorithm to fit 
FPs of the form 

N z 

logRo = alogcro + 51og/ 0 + c + , (1) 

i=l 

where for N z = 0 one recovers the standard form with free 
parameters a, b, c. Using N z = 0 and U max weights, defined 
as the inverse of the comoving volume in which a given 
galaxy would be observable, we initially follow the analy¬ 
sis of S13. However, we do not adopt the iterative clipping 
of S13, which is not justified as the distribution of residuals 
does not feature significant outliers or prominent tails. Con¬ 
sequently, the resulting root mean square residual of our fits 
is about 8 % larger (see Tabic [TJl , while there is good agree¬ 
ment in the best-fit coefficients with those in S13. Note that 
the analytic expressions for the errors on the fit parameters 
given in S13 are only applicable if measurement errors in cto 
and Jo were negligible compared to those of Ro, and if in¬ 
verse variance weights were used, neither of which holds true 
in their, and our, analysis. Instead, we estimate parameter 
errors from a repeated analysis of 100 bootstrap realisations 
of the underlying catalogue. 

Umax weights allow for the study of an ensemble of 
galaxies by removing the effects of Malmquist bias. Under 
the assumption that the FP and the galaxy selection process 
is redshift-independent except for the loss of faint galaxies at 
higher redshift due to the flux limit, the use of U max weights 
essentially fills in the loss of fainter galaxies at high redshift 
by proportionally upweighting them at low redshift. The FP 
derived in this way is a fair estimate of the FP for the en¬ 
semble under the aforementioned assumptions. 

However, for cosmological applications we require the 
FP to provide an unbiased estimate of galaxy size for each 
individual galaxy. We quantify the relative deviation of the 


2 http://sdss.physics.nyu.edu/vagc/ 

3 http://gax.shao.ac.cn/data/Group.html 
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Figure 2. The fundamental plane (FP) in the r-band used in 
this work (unit weights, including redshift dependence), shown 
for a random subset of the full S13 sample. The blue contour 
line encompasses 90 % of the galaxies in the most dense part of 
the distribution. The red line indicates a one-to-one relation. The 
inset shows the distribution of FP residuals A as defined in Eq.(2] 
which is well described by a Gaussian (red curve). 


galaxy size Ro from the expectation given by the FP, de¬ 
noted Rfp, via 


A = In 


Ro 


Ro 


Rfp Rfp 


-l + O 


( Ro 
\-Rfp 



( 2 ) 


which is an unbias ed estimator o f weak lensing convergence 
to first order (see iHuff fe Graved 1201 il l. We find A to be 
Gaussian distributed to good approximation (see Fig. H. 
The average of A in narrow redshift bins is shown in Fig. [3] 
While the mean of A over the whole sample vanishes by con¬ 
struction, the standard FP with U max -weighting and without 
explicit redshift dependence shows a strong scaling with z, 
clearly rendering it unusable for cosmological applications, 
particularly the extraction of peculiar velocities from ob¬ 
served redshifts. Additionally, our attempt to measure in¬ 
trinsic size correlations using the standard FP would neces¬ 
sarily show a large signal even in the absence of true intrinsic 
size correlations, just because correlation functions of galax¬ 
ies always pair up those that are nearby in redshift, meaning 
that the average residual for the galaxies in a pair will always 
have similar size and the same sign, and thus contribute a 
positive correlation. 

The redshift scaling of FP residuals is mainly caused by 
the flux limit of the SDSS Main sample which entails that 
we miss galaxies with small apparent size at higher redshift 
as they become too fainlQ. Moreover, some physically large 
galaxies at low redshift do not make it into the sample be¬ 
cause they become too bright or extended to allow for a 
reliable photometry or morphology measurement. As a con¬ 
sequence, we systematically miss galaxies above the FP at 


4 Note that there is also an implicit lower size cut in the S13 
sample, most likely driven by the requirement that they need to 
be reliably classified by Galaxy Zoo. However, we do not see any 
indications for a redshift trend in this selection effect. 
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Table 1. Fundamental plane (FP) fits according to Eq. [T] for the r and i bands. The fits differ by the choice of Vmax or unit weights, 
and by the inclusion of an explicit redshift dependence (‘zdep’). The rms deviation from the FP is given by opp. 


model 


band 


a 



b 



c 



di 

d2 


d<3 

crpp 

Vmax, no 2 dep 

r 

1.022 

± 

0.008 

-0.740 

± 

0.005 

-7.79 

± 

0.04 








0.100 




i 

1.050 

± 

0.008 

-0.745 

± 

0.005 

-7.78 

± 

0.04 








0.099 

no 

Vmax 

, no zdep 

r 

0.925 

± 

0.003 

-0.715 

± 

0.001 

-7.33 

± 

0.01 








0.084 




i 

0.950 

± 

0.003 

-0.727 

± 

0.001 

-7.37 

± 

0.01 








0.082 

no 

Vmax 

, zdep 

r 

0.784 

± 

0.003 

-0.682 

=t 

0.001 

-6.86 

± 

0.01 

1.75 

± 0.11 

-9.28 

± 

1.00 

28.6 

± 2.8 

0.077 




i 

0.816 

± 

0.003 

-0.696 

± 

0.001 

-6.93 

± 

0.01 

1.64 

± 0.11 

-9.50 

± 

1.00 

31.1 

± 2.9 

0.076 


low redshift and below the FP at high redshift, which leads 
to the tilt in the average FP residual A seen in Fig. [3] The 
Fmax weights exacerbate this trend because they upweight 
small low-redshift, and thus faint, galaxies. 

Note that a ph ysical reds hift evo lution o f FP parameters 
was observed bv lFernandez Lorenzo et alJ (1201 il l in a small 
sample covering the redshift range from 2 = 0.2 to 2 = 1.2. 
Within the narrow redshift range of the S13 sample, such 
trends ar e dwarfed by th e select ion effects discussed above. 
However, iBernardi et al.l (1200381 1 claimed to have observed 
weak redshift evolution in their FP, on top of selection ef¬ 
fects caused by the apparent magnitude limit, derived from 
a small SDSS sample over a redshift range compatible to 
ours. 

Switching to unit weights indeed reduces the redshift 
scaling of the mean residual (see Fig. 0. but only by incor¬ 
porating an explicit redshift dependence with N z = 3 in the 
FP relation of Eq. 0 does the average A become consistent 
with zero over the whole redshift range. As a welcome by¬ 
product, the rms deviation from the FP is now 18 % less 
than for the S13 results, without any manipulation of the 
residuals via clipping (see Table 0. A qualitatively simi¬ 
lar im provement was obtained bv lFernandez Lorenzo et al.l 
d201lll . The best-fit values of a , b, and c change substantially 
and now correspond to the FP extrapolated to 2 = 0. Note 
that the statistical uncertainty on the three classic FP pa¬ 
rameters remains unchanged, which suggests that they are 
nearly uncorrelated with 2 as the fourth parameter in the 
fit. Figurc[2] shows our optimal r-band FP, with unit weights 
and including the explicit redshift dependence, which we are 
going to use in the subsequent analysis. 


4 TWO-POINT CORRELATIONS 

We are interested in spatial correlations of the FP residuals, 
A, with each other and with the density field. To interpret 
these measurements, we also require the galaxy bias in order 
to translate a correlation with galaxy positions into a cor¬ 
relation with the matter distribution. We employ a suite of 
two-point correlation functions, measuring galaxy clustering 
(£gg), auto-correlations of A (£aa), and a cross-correlation 
between galaxy positions and A (£gA), with the following 
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2 

Figure 3. Mean relative deviation from the FP, A, as a func¬ 
tion of redshift. Blue diamonds correspond to a FP derived with 
Vmax weights and no explicit redshift dependence, as in S13. Red 
dots (black squares) show A for a FP with unit weights and with¬ 
out (with third-order) redshift dependence. Error bars are shown 
throughout but may be smaller than the symbols. 
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estimators 


£gg = 
CgA = 
£aa = 


DD — 2 DRu + Rr>Rr> 
RnRn 

F\D — F\Ru 
RfRo 
F xF x 
RfRf 


(3) 

(4) 

(5) 


We have suppressed the arguments of the correlation func¬ 
tions which are calculated over a range of (r p , n) bins, where 
r p is the comoving transverse galaxy pair separation and n 
is the comoving line- of-sight separatio n. The Landy-Szalay 
clustering estimator dLandv fe Szalavlll993l ) uses combina¬ 
tions of galaxy pair counts, denoted by XY for two samples 
X and Y. Here, D corresponds to the density tracer samples, 
i.e. either S13 or VAGC, while Rx denotes sets of random 
points with the same selection function as sample X. The 
terms involving FP residuals are given by 


FxX = 

Ai 0'l*>; 

(6) 


ieFjex 


FxFx = 

J2 AiA t0i*>, 

(7) 


ieF.jeF 


where (j\i) is a selection function that is unity if galaxies i 
and j lie within the respective (r p , n) bin and zero otherwise. 
Note that the subtraction of the FP residual count around 
random points in Eq. @ suppresses any additive system- 


© 2015 RAS, MNRAS 000. IRQ 



























Detection of correlations of fundamental plane residuals 5 


atic contribution to the correlation function £ g \. To reduce 
the shot noise contribution from randoms, we use random 
catalogues which are ten times larger than the real dataset 
when they correspond to the S13 sample, and twice larger 
than the real dataset when they correspond to the VAGC 
sample. The resulting pair counts need to be re-normalised 
accordingly, which we have not explicitly shown in Eqs. m 
to (0. 

We calculate the projected correlation functions via in¬ 
tegration over II, 


/ Ilmax 

£ab(r p , n) dn , (8) 

-n max 

which in practice is approximated by summation over II bins 
of width An = 10h^ 1 Mpc, with n max = 100h _1 Mpc. We 
also compare the clustering signals of galaxies above (A > 0) 
and below (A < 0) the FP, using the statistic 


Agg ( r p) 


w gg( r pi A > 0) — i 
w gg ( r pi A < 0) 


(9) 


A significant deviation of A gg from zero indicates a system¬ 
atic variation of the clustering properties across the FP. 

To get the covariance matrix for the correlation function 
measurement, we divide the samples into 100 jackknife re¬ 
gions of approximately equal area on the sky. Each region is 
about 64 square degrees in area (~ 8 degrees on a side). The 
size of the jackknife regions constrains our ability to measure 
the covariance matrix at large physical scales. To overcome 
this problem, we impose a redshift cut, z > 0.05, on our 
sample and also limit ourselves to scales r p < 30h^ 1 Mpc 
when interpreting the signals. 

While jackknifing the sample, we discard approximately 
8% of the S13 sample due to bright star masks and elim¬ 
ination of survey edges. After masking and imposing the 
redshift cut, we are left with 80,550 galaxies in the S13 sam¬ 
ple and 423,053 galaxies in the VAGC sample that we use 
in the final measurements presented in the following sec¬ 
tion. While these masks are uncorrelated with large-scale 
structure, they could in principle degrade the completeness 
and purity of the environment classification by spuriously 
removin g whole gal a xy gr oups or group members from the 
sample. Yang et ah, |2007) correct for the impact of survey 
edges, and our masking actually alleviates this effect further 
by removing the survey regions with the most complex ge¬ 
ometries. Bright star masks only account for a few per cent 
of the total area masked, so that their impact is negligible. 

We find strong correlations between the correlation 
function measurements in different r p bins on all scales con¬ 
sidered, in particular for tn gg , so that the full jackknife co- 
variances are used throughout. When inverting these covari¬ 
ances for likelihood analysis, we account for the noise due to 
the finite number of jackknife samples by employing an ap- 
proximat ely unbiased est imator for the inverse covariance, 
following iKaufmanl d 19671 '). 

We carry out the correlation function measurements 
separately in the r and i bands. The two-point functions 
in the two bands are highly correlated, so there is not much 
extra cosmological information when adding another pass- 
band, but requiring consistency of the results is a useful 
sanity check. 


5 SIGNALS AND MODELLING 
5.1 Clustering 


We measure the galaxy clustering signal of our density tracer 
samples with the primary goal of obtaining an effective lin¬ 
ear and dete rministic galaxy bias , 6 g . Closely following the 
formalism of lBaldauf et al.l J20ld '). we construct the follow¬ 
ing model for the projected correlation function, 


s (r r )-2bl Will f 

Jo l=0 L °s J Jo 

x £ 55,21 {yjx 2 +r%,z S ) L 21 


max 

dx (10) 


Vx 2 + 1 


+ Cic , 


where we introduced the growth rate f(z) = 
d [lnD( 2 )] /d [In a], with D(z) the growth factor and 
a = 1/(1 + z) the cosmic scale factor. The L 21 denote 
Legendre polynomials, and the < 22 1 are polynomials in 
the variable /3 = f(z)/b e whose explicit form is given in 
Eqs. (48) to (50) of iBaldauf et all (120101 ). We relate the 
multipoles of the matter correlation function to the matter 
power spectrum, Pg, as follows, 


£s5,2i(r , z) = 


(-1)* 

27T 2 


/” 


dk k 2 j 21 (kr) Ps(k,z) 


(ii) 


where the j 2 i are spherical Bessel functions. Throughout, 
we employ the full non-linear matter power spectrum as a 
means to phenomenologically extend our models into the 
mildly no n-linear regime, using the updated halofit non¬ 
linear fit_JTakahashi^t_aL[j2013) and the transfer function 
fit bv lEisenstein fe Hu 1 19981 ). 


We average the correlation function model over redshift 
with the weight 


W(a) 


Pa(z) Pbjz) \ f Pa(z) Pbjz) ' 

\ 2 {z)x'{z) [J ~ X 2 (z)x'(z)_ 


( 12 ) 


as motivated and derived in iMandelbaum et al.l (1201 ll ). The 
prime denotes the derivative of comoving distance with re¬ 
spect to redshift. We have written the weight more generally 
for the cross-correlation of two samples a and b, with red¬ 
shift probabilit y distrib uti ons p a (z) andpb(z), respectively. 
We follow iRoche fc Ealesl (1 19991 ) in estimating the integral 
constraint, Cic, via the random pair counts. Since the pro¬ 
jected correlation function is dominated by contributions of 
pairs with small line-of-sight separation, we find it is a fair 
approximation to calculate Cic using the projected matter 
correlation function, W5s(r p ), in combination with the ran¬ 
dom pair counts stacked along the line of sight. Moreover, 
the contributions by multipoles beyond l = 0 in Eq. m are 
small and hence can safely be neglected in the calculation 
of the integral constraint. Consequently, we compute 


Cic 


J2 rp [RdRd] (r p ) W5s{r P ) 

Er p [-Rd-Rd] (r p ) 


(13) 


where W55 is given by the monopole of Eq. (m divided by 

bl 

Using random pair counts over the range r p = 
[0.2; 200] Mpc /h, we obtain Cic = 0.29 Mpc /h for the NYU 
VAGC sample and Cic = 0.26 Mpc /h for the S13 sam¬ 
ple. The values for Cic are at least an order of magnitude 
smaller than w gg at the largest scales we consider. We then 
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Figure 4. Top : relative difference between the clustering corre¬ 
lation of galaxies above the fundamental plane (A > 0) and below 
the fundamental plane (A < 0), as a function of comoving trans¬ 
verse galaxy separation r p . Results for the r-band (black dots) 
and 2 -band (grey circles) are very similar. Centre/bottom: galaxy 
density-size correlation function, w g \ (centre), and galaxy size- 
size correlation function, w\\ (bottom). The best-fit model for 
each individual signal is shown as the red solid line. The vertical 
lines delimit the r p range used for the fit. Note that the correla¬ 
tion functions have been rescaled by a factor r p . Throughout, the 
S13 sample was used as the density tracer. 

fit the model of Eq. (nnj with fag as a single free parame¬ 
ter, using four bins over the range r p = [5.6; 30] Mpc /h. 
This results in a reduced \ 2 close to unity, with best-fit pa¬ 
rameters 6 g = 1.18 ± 0.09 for the NYU VAGC sample and 
b g = 1.37 ± 0.09 for the slightly more luminous S13 sam¬ 
ple (erro rs are la), in good agreement with previous results 
(see e.g. Zchavi et alj|201l] who also worked on SDSS Main 
samples and used a similar fit range). 

5.2 FP residual size correlations 

In Fig. [4] we show results for the spatial correlations of FP 
residuals A using the statistics A gg , w g \, and w\\ as de¬ 
fined in Sect. [4] The corresponding significance of the sig¬ 
nals, each measured over the range r p = [0.2; 30] Mpc /h 
and determined using the full covariance matrix of the data 
points, is given in Table [2] We observe a significant negative 
correlation in both A gg and w g \, plus a tentative positive 
signal in w\\, with no dependence on the passband in which 
the FP was determined. Using the S13 sample also for the 
galaxy number densities, produces marginally stronger de- 


Table 2. Significance of correlation signals in multiples of the 
width a of a Gaussian distribution, measured over the range 
r p = [0.2; 30] Mpc /h. Columns indicate whether the S13 or 
VAGC samples were used as the density tracer, and whether 
the r or i bands were used to determine the fundamental plane. 
The bottom half of the table lists systematics tests, re-computing 
the statistic with reshuffled sky coordinates (‘rand, pos.’) or us¬ 
ing only galaxy pairs with line-of-sight separation in the range 
II S [200; 400] Mpc//i (‘large IT). 


signal 


significance ( xcr ) 



S13, r 

S13, i 

VAGC, r 

VAGC, i 

™gA 

3.80 

4.18 

3.17 

3.07 

wx\ 

2.71 

2.36 

2.71 

2.36 

Agg 

4.28 

4.40 

3.35 

3.92 

™gg. large n 

1.51 

1.25 

1.02 

1.02 

Wg\, large n 

0.34 

0.52 

0.82 

0.95 

large II 

0.67 

0.83 

0.67 

0.83 

w g \, rand. pos. 

0.09 

0.06 

0.43 

0.27 

W\\, rand. pos. 

0.23 

0.37 

0.23 

0.37 


tections because it traces slightly more massive structures 
on average (see Fig.[TJ). A negative w s \ correlation implies 
that galaxies below the FP, i.e. with negative A, are prefer¬ 
entially located in environments with higher density. This is 
in line with the negative clustering ratio, A gg , which means 
that the galaxies below the FP are more strongly clustered, 
again pointing towards a high-density environment for these 
objects. 

To further elucidate this environment dependence, we 
compute the mean FP residuals for the different environ¬ 
ment subsamples introduced in Sect .[5] As is evident from 
Fig. [5] the brightest galaxies in groups (BGGs) tend to lie 
well above the FP, while satellite galaxies have preferentially 
negative A, with both trends showing only mild redshift evo¬ 
lution over the range we can probe. Isolated (‘field’) galax¬ 
ies constitute by far the largest subsample (see Section [2]). 
Their mean A is driven towards slightly negative values, in 
particular at low redshift, to compensate for the strongly 
positive A of BGGs because, by construction, our FP forces 
the mean A of the full sample to zero in every redshift bin. 
Since field galaxies and satellites dominate the contributions 
to w g \ because of their numbers, the resulting correlation is 
negative. 

The trends in Fig. [5] can also explain the result of neg¬ 
ative Agg. Since a halo can by definition host at most one 
BGG but arbitrarily many satellites, the latter are expected 
to reside on average in more massive haloes which there¬ 
fore are upweighted in correlations of satellites. Indeed, the 
haloes of BGGs in the S13 sample have an average number 
of 4.4 members while for satellites this number is 16.7. Con¬ 
sequently, we expect the satellite subsample to have higher 
galaxy bias than the BGG subsample. Since they are of ap¬ 
proximately equal size, galaxies below the FP should cluster 
more strongly than those above the FP, leading to A gg < 0. 
Field galaxies do not drive the effect because they are more 
weakly clustered, especially on the mildly non-linear scales 
for which we detect a signal, and, besides, are more sym¬ 
metrically distributed above and below the FP. 

In summary, the spatial correlations of FP residuals 
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we observe appear to be linked t o a dependenc e of the 
FP on the galaxies’ environ ment. iBernardi et al.l ([2003a) 
and lLa Barbera et, all l|201ij i both investigated the depen¬ 
dence of FP parameters on estimates of local density, using 
large SDSS samples at low redshift. Their results can be re¬ 
interpreted as the FP residual, A, smoothly increasing with 
density, in line with our findings for BGGs. Analogous trends 
have also been ob served for the central galaxies of more mas¬ 
sive s ystems fe.g. IBernardi et al.ll20071 : lvon der Linden et alJ 
1200711 . 


However, none of these works considered satellite galax¬ 
ies separately. It is suggestive to relate the consistently neg¬ 
ative A for satellite galaxies, as seen in Fig. [5] to systematic 
trends in galax y size at fixed mass. Evidence for such tre nds 
is either weak jR.ettura et al.l l201fljj_ Maltbv_et al. 201Cll i or 
opposite to the effects we see (Papovich et ah 2012h . but 
note that these works used far smaller samples. Moreover, 
trends in A could equally well be caused by a density depen¬ 
dence of the velocity dispersion , 00 , or the surface bright¬ 
ness, To - iBernardi et al.l (I2003bl l found that, in high-density 
regions, 00 and Rp increase while I p decreases, althou gh all 
dependencies were weak (see also iGooper et al.l 120121 for a 
similar result on the Rp dependence on density as measured 
in DEEP2/3). 

Systematic deviations from the FP have also been inter¬ 
preted as being d ue to variations in the mass-to-light ratio, 
M/L, of galaxies dGraves fe FabeJl20ld h lLa Barbera et alJ 
(f2010j) accordingly concluded from their FP analysis that 


galaxies in high-density environments have higher M/L on 
average. We can confirm this trend for BGGs and addition¬ 
ally infer that satellite galaxies have lower M/L than corre¬ 
sponding galaxies in the field, in agreement with the picture 
of mass stripping of galactic subhaloes falling into larger 
group and cluster haloes. 

Note that we have repeated the FP fit of Section[3]with 
either a local galaxy number density estimate or the distance 
to the nearest neighbour galaxy as an additional variable, 
but have found neither a change in the best-fit FP parame¬ 
ters nor a decrease in the scatter of the FP residuals. This 
implies that these two variables are not sufficiently good in¬ 
dicators of the environment dependence of the FP that leads 
to the spatial correlations observed in this work. 

We propose a simplistic model for the spatial correla¬ 
tions of the FP residuals, assuming that A is proportional to 
the matter density contrast, S, 


X(x) := B S(x) . 


(14) 


There is currently no clear motivation for a physical 
ansatz; this c hoice is just one of the simplest possible. 
lAlsing et al.l ll2014h used an analogous approach for mod¬ 
elling g alaxy size res id uals not related to an FP (see also 
ICiarlariello, Crittenden fe Pace! 120141 for a more sophisti¬ 
cated ansatz). Using Eq. (1141) . the projected correlation func¬ 
tions read 


w s \(r p ) = b g Bw ss (r p ) ; (15) 

w\\(r p ) = B 2 wss{r p ) , (16) 


with 


0.15 


0.1 

~ 0.05 

0 

-0.05 

- 0.1 


Figure 5. Mean relative deviation from the fundamental plane, A, 
as a function of redshift for subsamples of S13 defined according to 
environment. The colour coding is the same as in Fig.[T] Note that 
this was obtained for the r-band fundamental plane fit without 
Umax weights and with redshift dependence included. 
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where Jp is the Bessel function of the first kind of order 
zero. Here, we have assumed that the chosen n max is suffi¬ 
ciently large that the line-of-sight integration in the model 
can be extended to infinity. Furthermore, we have neglected 
redshift-space distortions. Note that the redshift distribu¬ 
tions of the S13 and VAGC samples are sufficiently similar 
that we expect negligible differences between the effective 
galaxy bias as measured from the VAGC clustering signal 
and the galaxy bias entering the cross-correlation w s \ with 
the VAGC sample as the density tracer. 

We fit this model over the range r p = [2.4; 30] Mpc /h 
separately to the correlation functions with the amplitude 
B as the free parameter. Results are given in Table [3] and 
the best-fit model in each case is also over-plotted in Fig. [4] 
We refrain from including smaller scales as the assumption 
of an effective linear galaxy bias, as well as the model of 
Eq. (1141) . are likely to break down. This implies that the 
model parameter constraints are weaker than the original 
detection significance down to sub-megaparsec scales would 
suggest. 

The constraints for w s \ are independent of the passband 
and density tracer and yield B ~ —0.01 with about a 50% 
error (Ict; note the error on galaxy bias has been propagated 
but is negligible). The reduced y 2 is quite small in this case, 
which may indicate that the jackknife covariances provide 
very conservative errors in the r p range of the fit. In contrast, 
the model yields a bad fit to w\\ and prefers an amplitude B 
that is substantially larger (note this statistic is not sensitive 
to the sign of B). We speculate that this could be a sign of 
a failure of our model. Higher-order contributions from a 
non- linear relation between A and 5 or a density weighting 
(see iHirata fc~S eliak 2004 for an analogous calculation in 
the context of intrinsic galaxy alignments) could have more 
impact on w\\ than on w s \ and boost in particular the 
smallest scales included in the fit range. 


5.3 Systematics tests 


r°° r°° au 1 u | 

wss{r P ) = / dzW(z) / —-—— Jo{k±r p ) P s (k±,z) , 

Jp Jp zn 

(17) 


We also carry out several tests for systematic errors. The 
first test is to measure w g \ and w\\ with the line-of-sight 
integral in Eq. © ranging from 200 to 400Mpc/h instead, 
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Table 3. Constraints on the intrinsic size correlation amplitude 
B, using the range r p = [2.4; 30] Mpc//i. The constraints on w g x 
were obtained with both the S13 sample (S) and the VAGC sam¬ 
ple (V) as the density tracer. Note that wxx can only constrain 
the absolute value of B. 


signals 

r-band 

B (la) 

Xred 

z-band 

B (la) 

X?ed 

w gX (S) 

-0.012 ±0.005 

0.08 

-0.010 ±0.005 

0.02 

™gA (V) 

-0.009 ± 0.006 

0.18 

-0.009 ± 0.005 

0.08 

w xx 

±0.038 ± 0.006 

2.81 

±0.035 ± 0.006 

2.43 


considering both positive and negative values of II (referred 
to in plots and tables as the ‘large II’ test). On these large 
scales, true intrinsic FP residual correlations should be neg¬ 
ligible, so this measurement can reveal whether there may be 
spurious signals due to image processing systematics, such as 
misestimation of the sky level around bright galaxies (which 
affects the flux and size measurements of nearby galaxies 
out to 100"; see, e.g.. lAihara et ahll201ll for estimates of the 
severity of the effect in DR8 data). Note that for very large II 
we expect non-negligible measurements of w g x and wxx due 
to lensing magnification, but the separations used are small 
enough that we do not expect any detection of this effect. Ta¬ 
ble [2] lists the significances of correlation signals in the large 
II statistic, measured over the range r p = [0.2; 30] Mpc /h. 
The results for all two-point statistics are comfortably con¬ 
sistent with zero. Note that the measurements for different 
passbands and density tracers are highly correlated, so that 
similar outcomes for the detection significances are expected 
across the columns of Table [2] 

Another systematics test involves permuting the sky co¬ 
ordinates of the galaxies, i.e. randomising the (RA,DEC) 
pairs in relation to the redshift and FP residual. This test 
(referred to as ‘rand, pos.’) is particularly sensitive to a non¬ 
zero mean of A as a function of redshift; see the issues raised 
in Section [3] Neither w g \ nor w\\ show any signs of signif¬ 
icant signals in this statistic for our final unit-weight and 
redshift-dependent FP (see Table 0, although both were 
clearly non-zero for the ‘classic’ choice of FP discussed in 
Section[3] Note that we observe a low-level negative correla¬ 
tion of A with random points, which is constant as a function 
of r v . This systematic remains below the signal in absolute 
value on all scales and is subtracted off automatically by the 
estimator of Eq. Q. 

When computing the model for our correlation func¬ 
tions (see Eg. HT1). we assume II max —> oo in Eq. ©. This 
is justified since our choice of n max = 100/i -1 Mpc when 
computing the correlation function from the data is large 
enough to capture almost all of the information. Still, to 
test for the choice of II max and All, we repeat the anal¬ 
ysis with a) n max = 50/i -1 Mpc and All = 5h _1 Mpc, 
and b) II max = 200/i -1 Mpc and All = 20/U 1 Mpc. Using 
the r-band and the S13 sample as density tracer, we find 
B = -0.011 ± 0.005 for case a) and B = -0.011 ± 0.006 
for case b), fully consistent with the results for the default 
choice of II max and All (see Tabled. Hence, conclusions are 
not significantly affected by our particular choice of line-of- 
sight binning. 

It is known that the standard SDSS pipeline 


has issues with the estimation of sky background in 
crowded regions, and especial l y for large objects (e.g., 
lAdelman-McCarthv et all 120061 : lAihara et al.l l201li l. This 
could induce spurious environment-dependent galaxy size 
measurements a nd thus affect the spatia l correlations we 
wish to detect. iHvde fc Bernardl (l2009l 'l found that the 
SDSS pipeline parameters underestimate the total flux and 
size for large galaxies when compared to their own data 
reduction. However, note that these authors worked on 
DR6 while we use the galaxy parameters from DR8 for 
which t he treatment of larg e galaxies was substantially im¬ 
proved (lAihara et al.l 120111 1. None the less, we employ the 
magnitude and radius corrections in Eqs. (4) and (5) of 
IHvde fc Bernardl 120091 ) to re-fit an FP and re-measure the 
w g x and wxx statistics, as a conservative estimate of the 
influence of sky subtraction effects. We find that these cor¬ 
rections modify the correlation functions by much less than 
la on all scales considered, so our results are expected to 
be robust against such systematics. Finally, the consistency 
of measurements in the r and i bands also underlines that 
our analysis is not significantly affected by issues with any 
photometric measurements that affect these passbands dif¬ 
ferently. 


6 IMPACT ON COSMOLOGY 

The small scatter of galaxy radii around the fundamental 
plane has been exploited in several applications of large- 
scale structure cosmology. In the following we will estimate 
the bias incurred by two such applications if the spatial cor¬ 
relations between FP residuals detected in this work are not 
accounted for. 

6.1 Peculiar velocity power spectrum 

The peculiar velocities of galaxies probe the slopes of grav¬ 
itational potentials and thus can be used to constrain the 
growth rate of structure, f(z). If an independent distance 
indicator is available, the observed redshift, z, can be split 
into contributions by the smooth Hubble expansion, denoted 
by zh, and by the line-of-sigh t component of the peculiar ve¬ 
locity, u. The relation reads iHarriso 3!Hzi) 

(l + z) = (l + z H ) (1 + ^) . (18) 

The FP is a widely used distance indicator in peculiar ve¬ 
locity studies. It is assumed that any systematic offset in 
radius from the FP is caused by the peculiar velocity of 
the galaxy, where Rq is determined from the observed red¬ 
shift, and Rfp from zh- The translation from galaxy ra¬ 
dius to redshift is done via the angular diameter distance, 
Da (z) = y(z)/ (1 + z), w here \ denotes comoving distance; 
see I Johnson et ah] d2014l ') for details. 

The spatially correlated contribution to A, as defined in 
Equation ©. adds a systematic shift, Au, to the inferred 
line-of-sight component of the peculiar velocity, u 0 bs- We 
linearise the response of u 0 b s to A and write 

dli 

Mobs = u ± Au « u + — A = u±T(z)A, (19) 

A=0 

where T(z) encapsulates the response. Using Equation (1181) . 
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Figure 6. Predicted contamination by intrinsic fundamen¬ 
tal plane residual correlations of the peculiar velocity auto¬ 
correlation power spectrum, P uu (k,fi = 1), at z = 0.1. In the 
top panel the cosmological signal is given by the black solid line 
and the systematic contribution by the red dotted line. The best- 
fit value for B from the w g \ fit (r-band, S13 sample as density 
tracer) has been propagated. Its propagated la error results in 
the error envelope shown in grey. The bottom panel gives the ra¬ 
tio of the systematic over cosmological signal, with lcr uncertainty 
indicated by the grey dotted lines. 


and linking zh via Da(zh) and Rfp to A, we obtain the 
expression 


T(zu) = —c 


1 + z 
1 "b zh 


(1 + Zji ) 


x'(zh) 

X( z h) 


(1 + Zh) 


x'(~h) 

x(«h) 



( 20 ) 


where denotes the derivatives of comoving distance with 
re spect to redshif t . The first equality agrees with the result 
of I Johnson et all (I20l4 . Eqs. 7 — 9). Note that our expres¬ 
sion lacks a factor In 10 since we employ the natural loga¬ 
rithm in A rather than the decadic one, and that we find 
a different sign in the Jacobian of the mapping between A 
and lo g [Da{z)/Da{zh)\, in agreement with ISnringob et al.l 
(12014) . The approximation in the second equality is equiv¬ 
alent to assuming u/c 1, which is fair in the context of 
this forecast. 

Using linear theory, the peculiar velocity component 
u can be linked to the matter density contrast, 5 (e.g. 
iHamiltonl Il998l ) . An analogous link can be established for 
the shift A u by combining Equations I® and m, which 
results in 

u obs (k)=[~iH(z)(l + z)- 1 f(z)^ +T(z)Bp(fc), (21) 

where tildes mark Fourier transformed quantities, and fi is 
the cosine of the angle between th e wave vector k and the 
line-of-sight direction (see also iKoda et al.ll20l4 for a moti¬ 
vation of the first term in Eq. ED). Here, H(z) is the Hubble 
parameter. The ‘observed’ velocity power spectrum (con¬ 
taining the combined effects of peculiar velocity and intrinsic 
size correlations) then reads 

p:t\k, n,z) = | + t2 (z)B 2 } Pss(k,z ) . 

( 22 ) 

The first term corresponds to the standard velocity power 


spectrum (see e.g. iKoda et al.l[20l4) . while the second term 
marks a purely additive systematic contribution by intrinsic 
FP residual correlations. In Fig.|S]we plot the two terms and 
the relative strength of the contamination, for [i = 1 and in¬ 
cludin g the empirical damping term proposed bv lKoda et al.l 
d20l4) . We have propagated the best-fit value and error 
of B for the fit to ui g \ in the r-band and using the S13 
sample as the density tracer. Intrinsic FP residual correla¬ 
tions could constitute a significant systematic effect for the 
peculiar velocity power spectrum above the 10 % level for 
k > 0.04/i/Mpc and may dominate the signal according to 
our best-fit model beyond k = 0.1 h/ Mpc. While the statis¬ 
tical uncertainty on this statement is still large and a null 
signal for the contaminant is not rejected with high confi¬ 
dence (due to the limited model fitting range), this finding 
suggests that any cosmological analysis based on peculiar 
velocity measurements from the FP should at least test for 
the presence of this syste matic. 

I Johnson et alJ (120141 ) performed a cosmological peculiar 
velocity analysis using the 6dFGRS, measuring the veloc¬ 
ity power spectrum out to k = 0.15/i/Mpc. They found no 
indication for an excess signal on small scales, nor any dis¬ 
crepancies with the signals as derived from a sample of type 
la supernovae. However, the statistical error bars are suffi¬ 
ciently large to accommodate a systematic at the few tens 
of percent level. Note that, since the authors recast their 
measurements into apparent magnitude correlations, which 
is beneficial for statistical reasons, they may be prone to yet 
another systematic as it seems plausible that galaxies also 
exhibit correlations of their intrinsic fluxes. 

There are several avenues to mitigate any systematics 
arising from intrinsic size correlations. Apart from contrast¬ 
ing measurements with those using other distance indica¬ 
tors, one can check the dependence of the velocity power 
spectrum on /r. Since the systematic contribution derives 
from an effect in the plane of the sky, it does not depend on 
H and could thus be filtered out as a constant offset. More¬ 
over, cross-correlations between galaxy position and pecu¬ 
liar velocity should be helpful in calibrating out intrinsic 
size correlations. In this statistic the g alaxy-velocity power 
is purely imaginary (IKoda et al.l 120141 ). while intrinsic size 
correlations generate a real contribution, provided A can be 
modelled as a local function of the matter density contrast. 


6.2 Weak lensing magnification 

Probing weak gravitational lensing magnification via the 
spatial correlations of apparent galaxy sizes is regarded as a 
promising comple ment to the estab lished gravitational shear 
measurements (iHeavens, Alsing fc Jaffel 120131) . The signal- 
to-noise ratio of such measurements is limited by the in¬ 
trinsic scatter of galaxy sizes, which is substantially larger 
than f or t he corre sponding measure of galaxy ellipticity 
(I Alsing et al.ll20l4 ). Therefore it could be desirable to in¬ 
stead use the much smaller scatter of galaxy sizes around an 
FP, even if this could only be obtai ned f or a mu ch smaller 
early-type galaxy sample fsee lHuff fe Gravesll20l4 for a re¬ 
cent application). Intrinsic correlations of FP residuals con¬ 
stitute a systematic for this type of weak lensing measure¬ 
ment that is closely analogous to the issue of i ntrinsic galaxy 
align ments in gravitational shear statistics dTroxel fe Ishakl 

l20i4) . 


© 2015 RAS, MNRAS 000.1111121 

















































10 


B. Joachimi, S. Singh, R. Mandelbaum 


10 


-9 


1(T 


1.2 

1 

0.8 

0.6 

0.4 


100 


hamper a weak lensing shear measurement at these redshifts 
as well. It should be kept in mind though that, in order to 
make this prediction, we have used a simplistic model, ex¬ 
trapolated it substantially in redshift, and ignored a likely 
dependence on galaxy mass or luminosity. Clearly, however, 
intrinsic correlations of FP residuals need to be investigated 
in more detail before magnification measurements based on 
the FP become a viable cosmological probe. With a rea¬ 
sonable model for this systematic, weak lensing statistics in 
principle have the power to self-calib rate the different i ntrin- 
sic correlations, as demonstrated bv lAlsing et al.l (|2014l) . 


1000 


7 CONCLUSIONS 



Figure 7. Predicted contamination by intrinsic fundamental 
plane residual correlations of the tomographic convergence power 
spectrum, C KK (£), for the cross-power between two redshift bins 
of width 0.1 centred on z = 0.4 and 2 : = 0.8. In the top panel 
the cosmological signal is given by the black solid line and the 
absolute value of the systematic contribution by the red dotted 
line. The best-fit value for B from the w g \ fit (r-band, S13 sample 
as density tracer) has been propagated. Its propagated la error 
results in the error envelope shown in grey. The bottom panel 
gives the ratio of the systematic over cosmological signal, with la- 
uncertainty indicated by the grey dotted lines. 


We assess the impact of intrinsic size correlations on the 
convergence power spectrum, 

? w (x) q (j) {x) 


C^(£) = 


r 

Jo 


d\ 


X^ 


P 55 


X 


,x 


(23) 


which is a statistic that can be derived from estimators of 
the form given in Equation ©. The integration limit Xhor is 
the comoving distance to the horizon, and the lensing kernel, 
q, is given by 


a w (x) = 




2c 2 o(x) 


C Xhor ... 

/ dx s p W (Xs) 


4 

V Xs J 


(24) 

where p^\x) denotes the probability distribution of comov¬ 
ing distances for a galaxy sample i. We choose to investigate 
a cross-power between two redshift bins of width 0.1, cen¬ 
tred on z = 0.4 and z = 0.8, respectively. This signal is 
affected by correlations between the intrinsic sizes (or their 
FP size residuals) with the local matter distribution, which 
in turn contributes to the lensing effect on the background 
galaxies. This matches closely what we have measured via 
the ui g A statist ic. T he corre sponding systematic signal then 
reads (see also lAlsing et al.ll2014l ) 


cmo.Br-, x ^MPM p,,a, x 

Jo X" VX 


(25) 


where Equation m was used once more. 

Employing again the best-fit value and errors from the 
r-band fit to w g \, using the S13 sample as density tracer, 
we obtain the power spectra shown in Fig. [7] On all angu¬ 
lar frequency scales that are accessible to measurement and 
modelling we find potentially severe contamination of the 
same order as the cosmological signal. Invoking the anal¬ 
ogy with intrinsic alignments, this is perhaps not surprising 
because bright elliptical galaxies as found in the S13 sam¬ 
ple also show the strongest alignment signals which would 


Using a large sam ple of elliptical galaxies (defined as in 
ISaulder et alJ l2013l ) in the SDSS Main galaxy sample, we 
report a 3 — 4<r detection of the spatial correlation function 
of fundamental plane (FP) size residuals with the large-scale 
density field. This detection is robust to the passband used 
to make the flux and size measurements (r or i) and to the 
choice of density tracers (the elliptical galaxy sample or the 
entire Main galaxy sample). A null detection of several sys- 
tematics signals, such as correlations at large separations 
along the line-of-sight or with shuffled sky positions for the 
galaxies, rules out measurement-related systematics for our 
detection. 

For further insight into our findings, we note that when 
splitting up the galaxies based on their membership in 
groups, there is a clear trend for group central galaxies to 
have a positive residual (meaning their observed size ex¬ 
ceeds what is expected from the FP fits), and for satellites 
to have a consistently negative size residual. We argue that 
these dependencies are likely to drive the negative correla¬ 
tion between matter density and FP residuals that we ob¬ 
serve. Consistent with these trends, the clustering two-point 
correlation function of galaxies with positive FP residuals is 
lower than that for those with negative FP residuals, at ~ 4 a 
significance. 

Our findings can be interpreted as being caused by an 
above-average mass-to-light ratio for central galaxies and 
a below-average mass-to-light ratio for satellites, which is 
consistent with the standard picture for tidal stripping of 
satellites. In the group regime that we are mostly sensitive 
to there is multiple evidence in simulat ions for th e str ip¬ 
ping of dark matter haloes fe.g. lKravtsov. Gnedin fe Klypinl 
12004) as well as the hot gaseous component of galaxies 
( Kawata fc Mulchaevl 120081 ) during the infall into larger 
haloes. Observational results also point towards the exis- 
te nce of this effect; s ee for instance the direct measurements 
bv lSuvu fc Halkola i‘201 P : using strong lensing and Sifon et 
al. (in prep.) using weak lensing. Clearly, a more quantita¬ 
tive comparison between the effects of tidal stripping and 
the spatial correlations in the FP residuals is desirable, but 
we leave this to forthcoming work. 

We em phasise t hat the classic three-parameter FP as 
derived bv ISaulder et al. (i2013 ) may be useful to study the 
global properties of the galaxy sample, but cannot directly 
be employed in any cosmological application, nor in the work 
done in this paper. One reason is that the sample definition 
with 1 /Umax weighting is meant to remove the impact of the 
flux limit by upweighting the lower luminosity galaxies that 
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appear in the sample at low redshift but not higher redshift, 
so that the overall sample luminosity distribution and de¬ 
rived FP is not subject to Malmquist bias. However, when 
calculating correlation functions of FP residuals, or quan¬ 
tities derived therefrom, one always pairs up galaxies that 
are at nearly the same redshift, and thus it is critical that 
the residuals be zero at all redshifts, not just on average, in 
order to avoid finding a spurious positive correlation func¬ 
tion. The Vjnax weights conspire with an implicit lower limit 
on galaxy size, induced by the flux limit, to create a strong 
trend of increasing FP residuals as a function of redshift. 
We correct for these selection effects by explicitly including 
redshift as a fourth parameter in the FP fit. 

Our detection has major implications for two types 
of cosmological measurements. First, lensing magnification 
studies that rely on galaxy size residuals as measured from 
the FP may be affected at the tens of percent level. Within 
our relatively large error bars we even cannot rule out 
that the intrinsic correlation attains the same amplitude 
as the cosmologic al weak lensi ng signa l o ver a wide range 
of scales. ICiarlariello. Crittenden fe Pacel (l2014Tj explored a 
halo model in which galaxy intrinsic size correlates with the 
density Held, which, however, is not directly comparable to 
our measurements because velocity dispersion and surface 
brightness variations also contribute to the FP residual cor¬ 
relations. Interestingly, they reached similar conclusions for 
the level of contamination expected for the size-convergence 
power spectrum as we did in Fig. [7] 

It will be elucidating to explore if, and to what extent, 
the correlations we detect can be attributed to effects of 
galaxy size, and to attempt a direct measurement of galaxy 
size intrinsic correlations, although the much increased in¬ 
trinsic size noise contribution will make this challenging. Al¬ 
though it may appear that the predicted levels of astrophys- 
ical systematic could hinder_the c osmo logical exploitation of 
size magnification, lAIsing et alJ (120141 1 showed that a joint 
analysis of shear and magnification of galaxy sizes can lead 
to substantial improvements in the constraints on dark en¬ 
ergy, of order 50 %, even when simultaneously marginalising 
over flexible models for intrinsic galaxy alignments and size 
correlations. Stu dies analo gous to galaxy-galaxy lensing, as 
conducted bv lHuff fe Graved J2014l l. are not affected by this 
systematic as long as the lens and sourc e galaxy samples 
are well separated fsee lBlazek et al1l2012l for the equivalent 
study using galaxy ellipticities). Besides, the level of intrin¬ 
sic correlations of residuals in the photometric analogue of 
the FP, as employed in that work, is currently unknown. 

Second, attempts to use FP residuals as a distance indi¬ 
cator for the purpose of peculiar velocity studies could expe¬ 
rience systematics due to intrinsic correlations at the 10% 
level or above for wavenumbers k > 0.04 h /Mpc in the pecu¬ 
liar velocity auto-power spectrum. Our 3-4<r detection does 
not give strong constraining power on the amplitude of the 
effect, and our ad-hoc model choice may induce further un¬ 
certainty in the forecast. Therefore, both observational and 
theoretical effort is warranted to get a better understand¬ 
ing of spatial FP residual correlations and their impact on 
peculiar velocity cosmography. Fortunately, one can identify 
several mitigation approaches for this systematic, including 
the dependence on the line-of-sight angle p (since the in¬ 
trinsic correlations are in projection on the sky), the use of 
self-calibration by incorporating cross-correlations between 


the galaxy distribution and peculiar velocities, as well as the 
comparison with samples relying on other distance indica¬ 
tors. 

These findings underline that, despite past claims to the 
contrary, cosmology using galaxy surveys with today’s accu¬ 
racy requirements is impossible without a thorough under¬ 
standing of the underlying galaxy samples and their phys¬ 
ical characteristics. The simultaneous inference on cosmo¬ 
logical parameters and astrophysical effects related to in¬ 
trinsic galaxy properties, such as the correlations found in 
this work, is therefore likely to be the default mode for the 
analysis of the forthcoming large spectroscopic and imaging 
surveys. 
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